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I We re-express the Kerr metric in standard Bondi-Saches' coordinate near null in- 
finity . Using the uniqueness result of characteristic initial value problem, we prove 

^ . the Kerr metric is the only asymptotic flat, stationary, axial symmetric, Type-D solu- 

qh| tion of vacuum Einstein equation. The Taylor series of Kerr space-time is expressed in 

5^ ' terms of B-S coordinates and the N-P constants have been calculated. 

keywords: Kerr solution. Uniqueness Theorem 

> ■ PACS code : 04.20.-q, 04.20.Ex, 04.20.Ha 
^ . 

^ ■ 1 Introduction 

T— I ■ 

. After the work by Bondi et.al.[l], it is well-known the Bondi coordinates is a very natural 
^ ' choice when we want to describe the asymptotic behavior of gravitational field near null 
infinity X"*". Based on works by Penrose, Newman and Unti[2, 3], there is an elegant way to 
^ ■ re-express Bondi's work in N-P formulism. This also gives us a general formulism to describe 
the asymptotic structure of general asymptotic fiat space-time which is smooth enough near 
X"*". Using characteristic initial value (CIV) problem method, many authors[8, 9, 10] have 
shown the existence of null infinity in general case and pointed out the degree of freedom 
of gravitational field near null infinity. The CIV method has many advantages in dealing 
with gravitational radiation problem. Recently, this method has been used in numerical 
relativity. Winicour and his colleagues develop the CCM method. They want to combine 
the CIV method and standard Hamiltonian evolution method[ll]. On the other hand, Kerr 
solution is a very important exact solution of Einstein equation both in theoretical area 
and in application. An interesting question is whether Kerr metric describes the space-time 
outside a stationary rotating star. For long time, the Bondi coordinates of Kerr space-time 
is not very clear. For example, how to describe the Kerr space-time in Unti-Newman's 
general formalism[3]? The uniqueness theorem[12, 13] tells us that Kerr solution is the 
only asymptotic fiat, stationary, axial symmetric solution of vacuum Einstein equation with 
regular event horizon. From the application point of view, it is very difficult to get the detail 
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information about the event horizon of a space-time because of the infinite red-shift near 
horizon. An interesting question is how to identify the Kerr solution based on information 
near null infinity. This is more practical in future gravitational experiments. This idea 
can also be understood from the Geroch conjecture[14]. Obviously, stationary and axial 
symmetric condition is not enough because there are many asymptotic fiat exact solutions 
of Ernst equation. In next section, it is found such uncertainty comes from the homogeneous 
part of a control equation which comes from the Killing equation. The general solutions 
of that equation contain some free constants. These unknown constants are close related 
with Geroch-Hansen multi-pole momcnts[17]. In order to identify the Kerr solution, we use 
Petrov classification [4] and show that condition will help us to pick out Kerr solution finally. 
Further more, N-P constants for Kerr space-time is also calculated as a byproduct. 

This paper is organized as following : In section II, we prove an local uniqueness theorem 
of Kerr solution based on the information near null infinity. This theorem also tells us the 
standard Unti-Newman expansion of Kerr metric. The detail expression of this extension is 
contained in Appendix A. Appendix B contains some spin-weight harmonics which is useful 
for our calculation. 

2 Main theorem 

Let (M, g) be an asymptotic fiat space-time, (u, r, 0, (p) be the standard Bondi-Sachs coor- 
dinates. 

Theorem 1 Suppose {M,g) be an asymptotic flat, stationary, axial symmetric, Type-D, 
vacuum space-time with in a neighborhood of null infinity, then it is isometric to Kerr space- 
time in the Bondi coordinates neighborhood. 

In order to make the proof clearly enough, we divide it into two subsections. In first 
subsection, we calculate the Taylor series of general stationary axially symmetric space-time 
in Bondi coordinates. We find all Taylor coefficients can be expressed in terms of {"^q} 
and their derivatives. Unknown functions {^q} satisfy a linear inhomogeneous second order 
partial differential equation. We find the general solution of such equation has the form 
^0 = ^0 + -^'^ 2^fc+2,o, where is the special solution which corresponds to Kerr metric 
and D'' is a free constant, i.e. all axially symmetric and stationary solutions are characterized 
by the set of constant {D''}. This set of constants are close related with the famous Geroch- 
Hansen multipole. In section 2.2, with the help of Type-D condition, we show all {D'^} 
vanish and Kerr metric is the unique stationary, axial-symmetric, asymptotic fiat Tpye-D 
metric . 

2.1 Taylor series of general axial symmetric vacuum stationary 
spacetime 

Suppose (M, g) be a vacuum stationary axial symmetric space-time, t"- and 0" are two 
commutative Killing vectors. Near wc use the standard B-S coordinates to do the 
standard asymptotic expansion. The detailed construction of these coordinates is well-known 
and can be found in Re. [2, 3]. With this choice of coordinates, we also can choose a set of 
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null tetrad {l"" , n"" , m"" , m"'} , such that = (^-^^ and these tetrad are parallel-transported 
along Under such choice of coordinates, 0" = • The time-like Killing vector can 

be expressed in terms of null tetrad as t"- = Tl"- + Rn"" + Am"' + 74m". [t"", 0"] = means 
that T, R, A are independent of if. The null tetrad components of Killing equation for are 

-DR +{e + e)R + kA + kA = Q, 

-DT -{e + e)T -nA-nA- D'R + (7 + 7)/? + + = 0, 
-kT + TiR + DA + {e- e)A -6R+{a + p)R + pA + aA^O, 
-D'T - (7 + 7)r - uA- uA = 0, 

-tT + vR + D'A + (7 - 7)A -5T-{a + I3)T - /xA - = 0, 
-aT + \R + 5A+{a- (3) A = 0, 

-pT + pR + 6A- {a- P)A - pT + pR + 6A - {a - /3)A = 0. 

Using the Bondi Gauge k — e — ir — 0, p — p, t — a + P,we have 

-DR^O, (1) 

-DT - D'R + (7 + 7)i? + fA + rA^O, (2) 

DA - 5R + rR + pA + aA^ 0, (3) 

-D'T - (7 + 7)T - - zvA = 0, (4) 

-tT + PR + D'A + (7 - 7)A -5T -tT - pA-XA^O, (5) 

-(7T + Ai? + M+(a-/3)A = 0, (6) 

-pT + pR + SA- {a- P)A -pT + flR + 5A-{a- p)A = 0. (7) 

Here we use the standard notation of [2, 3, 4]. Differential operators in above equations are 
defined as 



D 



d_ 

/ d d d-d 

du dr d( 9^' 

It is well known that stationary solutions to Einstein's vacuum field equations are 
analytic [5]. Moreover, it is also known that asymptotically flat stationary vacuum solu- 
tions are not only analytic, but even admit an analytic conformal extension through null 
infinity[6, 7]. Keeping this result in mind, all geometric quantities (the coordinate compo- 
nents of null tetrad, N-P coefficients, components of time-like Killing vector and components 
of Weyl curvature) can be expressed in terms of power series of ^, for example 

T = T° + — + •••, 

r 

A = A" + — + ••• , (9) 
r 

some lower order Taylor coefficients of components of tetrad, N-P coefficients and components 
of Weyl tensor (up to 3th order) can be found in section 9.8 of [2]. 
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First of all, let's consider the function R. Eq.(l) and axial symmetric condition tell us 
that R = R{u,9). With the formal expansion of null tetrad and N-P coefficients [2, 3], zero 
order of Eq. (2) gives 

so R = R{0, (p). In order to get more information about R, higher order of Killing equation 
are needed. The first order of Eq.(3), Eq.(5) and Eq.(7) are 

SoR + A"" = 0, (11) 
+ - 5or° + - = 0, (12) 

2T°-R = 0, (13) 

where 5q = ^^^^ and " ■ " means (9„ . Because the space-time is stationary, there is no 

Bondi flux. This implies cr° = [2], then N-P equations tell us this leads ^'g = 0. Combining 
this condition with (11), (12) and (13), we get 

ii + A° = 0. (14) 

The second order of Eq.(3) is 

It has been shown that the right hand side of above equation is independent on m, so = 0. 
This implies A^ = and i? is a constant. So we can take R = 1 without lost of generality. 

With R = 1, Killing equations become 

-Dr+ (7 + 7) +fA + TA = 0, (16) 

DA + T + pA + aA^O, (17) 

-D'T - (7 + 7)T -uA-uA = 0, (18) 

-tT + P + D'A + (7 - 7)A - 5r - tT - //A - AA = 0, (19) 

-aT + \ + 5A+{a- (3)A = Q, (20) 

-pT + + M - (a - /?)A - pT + /2 + M - (a - /?)A = 0. (21) 

With Eq.(9), the non-trivial zero order Killing equation is 

T° = 0. (22) 

Non-trivial first order Killing equations are 

= 0, (23) 
-f^ + + = 0, (24) 

+ A^ + - = 0, (25) 

a" + 5oA° + 2a°A° = 0, (26) 
2T° - 1 = 0, (27) 
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which impUes = and = 0. Same order N-P equations also give ^3 = and ^4 = 0. 
Second order KiUing equations are 

Ti-l(^o^^rO) =0, (28) 

-2A^ = 0, (29) 
= 0, (30) 

^^^l + A'-SoT'^O, (31) 

lcT° = 0, (32) 

2T' - - a'a'' - ^0 - = 0, (33) 

where ^ is the spin-weight operator[2] and is defined as := {Sq + 2sa^)f, = — ^| ■ 

From these equations, we know ct° = 0, = 0, = 0, = i(^^ + '1'^), = = 0, 

= -|^\&2 + |'^o(^2 + "^'2)- It is worth to point out that the result cr° = tells us that the 
Bondi coordinates which is chosen as [3] are associated with the "good cut" of stationary 
spacetime, i.e. (T° = 0, so the freedom of super-translation has been removed. 

The third order killing equations are 

2T2 + i(^*5+ ^^0) = 0, (34) 

= 0, (36) 

SqA^ + 2a°A'^ = 0, (38) 

+ 1 + 1 ^^ro + SqA'' - 2a'' A^ + doA^ - 2a'' A" = 0. (39) 

Eq.(35),(38) imply 

= 0. (40) 

The spin- weight of ^'^ is 1, so it is a linear combination of spin- weight harmonics 

The axial symmetric condition implies m — 0. The behavior of spin-weight harmonic under 

the action of operators ^ and ^ are [2] 



P Im — V 2 s-lllm: 

oYlm — Ylm- (41) 

So we get \E'5 = c(u)iFi_o = c{u) sin 9. Detailed calculation on same order N-P equations also 
give 

= - ^3 = 0, = ^1 = ^1 = ^1 = 0, = (42) 
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Eq.(34) and = gives 

c^sin^ + c^sm6' = 2(c+c)cos6' = 0, (43) 

which imphes c + c = 0. The first order Bianchi identities tell us ^5 ~ ~ 0' which 
imphes ^2 — — ccos^ + C2 but it is well-known that ^2 — *2 f*^^ stationary case [2, 3], so 
c — 0. This means = Cisin^ and ^2 — ^2- The Komar integral shows — C2 is just the 
Bondi mass of the space-time. Additionally, Eq.(37) tells us that — 0. 

Forth order KiUing equations are 

3T3+(y + f ) =0, (44) 
4A=^=i^< (45) 

f' + ^{M+ ^*?)(*° + ^°) = 0, (46) 
i^oyi -l^-q,o + p^ + A^ + (*° + ^°)^2 ^ 3^3 _ ^^^3 ^ ^0^2 ^ 1^3 _ ^47^ 

+ + = 0, (48) 

+ ^4 ^ ^ + = 0. (49) 

where 7^ = -^^(aO^vl/O _ ^o^^O) _ i_f^o^ ^4 ^ _^^o^ ^4 ^ _|^2 ^ -i^^^o^ ^4 ^ 

^(^^[] + j^^*g), *i = -|*?*2 - If'^o (These results are got from same order N-P 
equations). 

The spin- weight of ^0 is 2, so Eq.(45),(48) imply 

*° = D^u) sin^ e, (50) 
Eq.(37),(45) eliminate time dependence of D^{u), i.e. = -D^sin^6'. Eq.(47) is 

= (51) 

Fifth order Killing equations are 

4r^ + (7^ + f ) - -^[A" - = 0, (52) 

= (53) 
o 

-2p' + 2r^ + (//^ + fl^) + ^ |^o^2 ^ ^^4 ^ ^^4 ^ (55) 

where 

p' = 0, 
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Eq.(52), (46) and Bianchi identities imply = zCisin^, Ci G R, where Ci is the Komar 
angular momentum. Eq.(53),(54) give 

^^ml + 5^1 = 10(*?)^ - 15*^*^. (56) 

The homogeneous part of above equation is 

^^*J + 5*J = 0. (57) 

Because spin-weight of is 2, it is a linear combination of {2^,0}- Using Eq.(41), the 
homogeneous equation is 

(_/2 _ / + 12) 2l^,o = 0, (58) 
which gives 1 — 3. The general solution of Eq.(54) is 

*° = (y - ^ + 2^,0- (59) 

(Bianchi identities insures ml = 0.) Bydefinition[2], the non-zero N-P constant for stationary 
axial symmetric space-time is 

Go^^iCif-dC^D'. (60) 

Until now, we have got series expression of tetrad components up to 4th order, N-P 
coefficients up to 5th order and Weyl components up to 6th order. To prove this theorem, 
all Taylor coefficients of all geometric quantities are needed. We use inductive method to 
solve this problem order by order. 

Suppose we have known Taylor coefficients of tetrad components up to (A; — 3)*^^ order, 
Taylor coefficients of connections up to (A; — 2)*'* order and Taylor coefficients of Weyl cur- 
vature components up to {k — 1)*'* order. The {k — 1)*'* order of KiUing equation (17) and 
(20) are 

-{k - 1)A''-^ + t''-' ^ ■ ■ ■ , (61) 
• • • + + ^A''-^ = 0. (62) 

where "• • •" means terms which only contain lower order coefficients. Based on the induction 
hypothesis, those terms are known. In order to solve these equations, we need coefficients 
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A*^ ^ and r*^ ^. From N-P equations, 



- 5*0 = -4q;*o + 4p*i = - ^ ^*oH . 

Dr = Tp + fa + *i ^ r'=-^ = -^^l^*J + ..., (63) 



Combining Eq.(61), (62) and (63), we get 



^^*J+(^i±ili^*J = .... (64, 



2 

The homogeneous part of above equation is 

+ (^±iM±l) ^ 0. (65) 

Because of Eq.(41) and axial symmetric condition, the general solution should be 

= + 2^+2,0 , (66) 

where is a special solution of eq.(64) and D'' is a constant. Obviously, Kerr solution 
satisfies all conditions of our theorem, so must exist. The concrete form of also 
can be got by direct calculation. One can express the "• • •" terms in Eq. (64) as a linear 
combination of spin-weight harmonics {2^,0}- The inductive method insures the maximal 
value of / in that expression will be finite for any given order, then we can get \E'o by 
comparing coefficients between bother sides of this equation. With the general solution of 
Eq.(63) will give t'^~^, a'^~^, X^~^ and ^'i. Further more, Cartan structure equations 
and Bianchi equations will help us to get other coefficients. 



Dp = p2 + |(7p 


=^ 


-{k- 




' ' ' ■) 


Da = ap + (3a 


=^ 


-ik~ 


2)a^-^ 


' ' 


DP ^ l3p + aa + *i 


=^ 


-{k- 






D'^2 - = 3p*2 - 2q;*i - A*o 


=^ 


-{k- 


3)*^ = 




D^s - 5*2 = 2p*3 - 2A*i 


=^ 


-{k- 


2)*^ = 


M + 


D*4 - 5*3 = p*4 + 2q;*3 - 3A*2 


=^ 


-{k- 






D^^Ta + fP + *2 


=^ 


-{k- 






Dp = pp + Act + ^2 


=^ 


-{k- 




= ^p'-' + ^2 + ---, 


Du ^ rX + fp + '^s 


=^ 


-{k- 


ly-' ■ 






=^ 


-{k- 


mi-2 -- 


~ ' ' ' ) 


D^ = pe + 


=^ 


-{k- 


3)et2 -- 


~ ' ' ' ) 


Dcu — pco + aco — {a + P) 


=^ 


-{k- 


3)a;^-2 


= -a^-'^ - P''-^ + ■■■, 


DX={a + P)e + {a + 0)^'' 


^ 


-{k- 


2)X'=-2 




DU = (a + (3)od + {a + P)uj - 7 - 7 




-{k- 


2)^fc-2 


— _7'=-i _ 7*^-1 + • • • . 



(67) 
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From above results, we find we can express all [k — 2)*^^ order coefficients of tetrad com- 
ponents, {k — 1)*^ order coefficients of connection components and /c*^ order coefficients of 
Weyl curvature in terms of ^'q, ^ derivatives of and lower order coefficients which we 
have known. The form of is given in Eq.(66). This means we can get all those coefficients 
for any given order. 



2.2 Uniqueness of Kerr solution 



In above subsection, we have got Taylor series of a general stationary axial symmetric metric. 
Prom Eq.(66), we can see that the freedom in each order Taylor coefficients are just the 
constant D'^ { k > 5). These arbitrary constants should be closely related to the famous 
Geroch-Hansen multi-pole moments[14, 15, 16, 17]. What we want to do in this section is to 
pick out the Kerr solution from those possible solutions, i.e. we need to fix value of {-D*^}- 
In order to do that, we consider the Petrov classification [4]. It is well known that the Kerr 
solution belongs to Type-D class, i.e. its Weyl curvature satisfies[4] 



K = ^1(^4)^ - 3^4^3^2 + 2(^3 



0. 



Write down the 12*'' order coefficient of above equation, we get 



(68) 



(69) 



In previous section, we have got 



*2 = C2, 



2^,0, 



*3 



2 
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(70) 



This constant is fixed in following way : from Komar integrals M — ^ Js^ *dt and J = 

Ma = ^ Js^ *d(f), we know = -M, *° = SiMa^ lYi^. Submit these into Eq.(70) 
then get 



*3 



3iMa 




-1^^2,0 



2^,0 



(71) 



Submit above result into Eq.(69), we find 



= 3M 



3iMa 
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3A/ 



3M^a' 167r 



2^2,0 



Solving the simple linear algebraic equation 3M 
can fix the value of is 



3M2a2 



'IGtt 



+ zV67rM^a. 



167r 
5 



(72) 



0, we 



(73) 



Submit above result into Eq.(60), it is easy to check the N-P constant of Kerr space-time is 
zero, which has been got by [19, 21]. 

In order to fix general D*^, wc also use inductive method and fix them order by order. 
Suppose we have known D'' up to order n. To get the value of we consider the (n + S)*'* 

coefficient of Eq.(68), a long but direct calculation shows it should be 



(74) 



here "• • •" also mean terms which only contain lower order coefficients. Prom Eq. (63), (66), (67), 
we know 



vn+l 



,T,n+l I ptn+l y 



n-3) 



n-2 
1 

n — 1 

1 -„ 



n+l 



+ 



+ 



^^^+1 + 



n 



(75) 



where \I'o^^ is the special solution of Eq.(64) which corresponds to Kerr solution. We have 
known Kerr solution belongs to Type-D, i.e. Eq.(74) holds for ^q, so Eq.(74) can be written 
as 



-3*1^0 



(n- 3)(n- 2)(n- l)n 



^^^"+1 + . . . = 0. 



(76) 



If the general ^q"*"^ in Eq.(75) also satisfies Eq.(74), i.e. 



-3*1*0 



(n-3)(n-2)(n- l)n 



-3*1*0 



(n-3)(n-2)(n- l)n 



r 2i^n 



+3,0 



+ 



0. (77) 



Because terms in "• • •" only contain lower order coefficients, they remain unchanged when we 
change to the general *r\ Obviously, ^—^^^^^^^^-^^ ^Y^ 

+3^0 is a non-zero function 

for any n, so the general solution of *o'^^ satisfies Eq.(74) means D"+^ = and *o''"^ = *o"^^ 
This tells us that Kerr solution is the only solution which satisfies all requirements of our 
theorem. 
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Remark : in above subsection, we proved our theorem under the requirement that the 
whole space-time is type-D. In subsection 2.2, we have shown that the freedom of vac- 
uum, stationary, axial-symmetric space-time are just a set of constants {-D^}- The rea- 
son why we need the condition type-D is to fix the value of {-D*^}- If the type-D condi- 
tion holds at several points in Bondi neighborhood and those points are not zero-points of 
(n-3)(n-2Xn-i)» 2^n+3,0) then, it is casy to see that {-D'^} should be zero, i.e. type-D condi- 
tion holds at several points will imply this condition holds in the whole neighborhood. This 
feature may be a practical method from the experiment prospective. That means only sev- 
eral points need to be checked for the type-D condition to see whether the space-time around 
us is Kerr space-time. This maybe a useful property for future gravitational experiments, 
such as "mapping space-time" [22] . 



Acknowledgement 

This work is supported by the Natural Science Foundation of China under Grant Nos. 10705048, 
10605006, 10731080 and K.C.Wong Education Foundation, Hong Kong. Authors would like 
to thank Prof. X.Zhang and Dr. J.A.Valiente-Kroon for their helpful discussion. 



Appendix A 



The asymptotic extension of Kerr space-time in B-S coordinates are 
1) Null tetrad 



r 



d_ 

dr ' 
d_ 

du 



-- + -~^(3cos^^-l) + 0(r-^^ 



d_ 
dr 



+ 



iMa 9 ^, 4, 
^cot- + 0(r-) 



a 



iMa 9 ^, A.] d 
■^oot-+0(r-) (78) 



3iMa . ^ Ma^ „^ a, 
sm9+ ' ^^---^^ 



d 



+0(r-)- + 



:i + CO 

V2r 



sin^ 9 + 0{r 

+ 0{r-* 



d_ 
dr 



d_ 



2) N-P coefficients 



a — 



a — 



- + 0(r-5), 
r 

3Ma^ sin 9 



2r4 

cot 61 3Ma^ sin 6* cos 6' 

+ ^^^^ + 0{r-^), 



2^/2r^ 



2V2r 

cot 9 SiMa sin 9 33 Ma'^ sin 9 cos 9 ^, 

+ + C>(r-^), 



2^/2r 2v^r3 



2^/2r 
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T = 

A = 

7 = 



SiMasinO ISMo^ sin ^ cos 6* 

+ + O r-^ , 



+ 0(r-^), 



2v^r 
Ma^ sin^ 6 

_1_ M SiMa cos e Ma^(3 cos^ ^ - 1) -5^ 
M {2y/2-l)3iMacose 3 Ma"^ {3 cos^ 9 - 1) -5^ 



2r2 

3iMa sin 6* Ma^ sin 6* cos 61 



4^4 



+ 0(r-^). 



4\/2r3 

3) Weyl curvature 

3Ma2sin2^ 



\/2r4 



^4 



+ 0(r- 



3iMosin^ 12Ma2 sin^cos^ 



+ Oir-% 



V2r5 



M 3iMacos^ 3Ma2(3 cos^ ^ - 1) ^, 



4r4 



3iMa sin ^ 



2A/2r4 



+ 



3z 6z 

^M^a sin 9 + —^Ma^ sin 9 cos 6* 

2^/2 ^2 



r"^ + 0(r" 



+ 0(r-^). 



(79) 



(80) 



Appendix B 

Some spin-weight harmonics 

1^1,0 



47r 



16% 



[cos9+l)e''f', 



, 3 . . 
— sm 9, 

Sn 



IGtt 
"3 
Stt 




(1 -cos^)e-^'^; 
sin^e*^, 



3 

' — cos^, 
An 



sin 9e 




167r 



(l-cos^)e^*, 
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Stt 



sin^, 



IGtt 



(l + cos^)e-*'^; 



3a 



967r 



247r 



167r 



247r 



(l + cos^)V^^ 
sin^(l + cos^)e^'^, 



sin^ e, 



sin^(l -cos^)e-''^, 



967r 



(1 -cos^)^e 



2 -2i<j). 



5 



247r 



sin^(l + cos^)e2^<^, 



247r 



(2 cos 6^ - 1)(1 + COS! 



3\l — sin^cos^, 



31 
3i 



247r 



247r 



IGtt 



(2cos^ + 1)(1 -cos^) 



sin ^(1 -008^)6-^^"^; 



sin^^e^^^ 



sin ^ cos ^e*"^, 

lOTT 



247r 



(3cos2^- 1), 



6^ 
3a 



167r 



sin 6 cos 



167r 



sin 6*6 



247r 



sin^(l -cos^)e^"^, 



247r 



(2cos^ + 1)(1 -cos^) 



— sin ^ cos ^, 
An 
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-i>2,_i = 3J^(2cos^-l)(l + cos^)e-^^ 



-2^2.2 = 3y'^(l-cos^)V^^ 



-2 



-2>2,-i = -3W^sin^(l + cos^)e-*^ 
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